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Метод оценки частоты генераторов в условиях непрогнозируемого изменения 
длительности интервала измерений* 


О. А. Сафарьян 


Исследуется метод оценивания частот одновременно и независимо функционирующих генераторов в условиях не- 
прогнозируемого изменения длительности интервала измерений. Рассматриваются погрешности, возникающие при 
использовании указанного метода из-за нестационарности частот генераторов на интервале оценивания частот. 
Первая составляющая ошибки связана с отклонением измеряемой фазы колебаний г енератора из-за собственной 
нестабильности частоты г енератора, вторая определяется непрогнозируемым изменением частоты генератора на 
интервале измерении, Отмечено, что уменьшение каждой из составляющих предъявляет взаимоисключающие тре- 
бования к длительности временного интервала. На основе известных соотношений, определяющих потенциально 
достижимое значение среднеквадратического отклонения частоты от номинального значения, получены выраже- 
ния, показывающие оптимальную длительность временного интервала измерении, В качестве критерия при выборе 
длительности временного интервала рассматривается минимум суммы двух ошибок. Приводятся аналитические 
соотношения, определяющие потенциально достижимые точности оценок частот. Доказываются несмещенность, 
эффект, ИвВНоСтТЬ И состоятельность получаемых оценок, 
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Введение. На современном этапе развития инфокоммуникационных технологий повышение точности 
и стабильности формирования частоты генераторов является одной из важнейших задач. Эти вопросы 
актуальны для систем связи радиолокации, радионавигации и метрологии. Однако во многих случаях, 
например, в радиолокационных или радионавигационных системах, можно ограничиться только зна- 
нием частоты генераторов, чтобы учесть это значение при определении навигационных параметров 
объекта на местности и в пространстве или при поверке приборов [1—5]. 

В настоящее время для определения и стабилизации частоты генератора наиболее широко ис- 
пользуется способ с применением фазовой автоматической подстройки частоты (ФАПЧ) [6-11]. Реали- 
зующее данный способ устройство включает эталонный и подстраиваемый генераторы, фазовый детек- 
тор и управляемый элемент. Однако высокостабильный генератор является сложным и дорогим техни- 
ческим устройством. Это связано с построением высокодобротных контуров и формированием высоко- 
стабильных колебаний [12, 13]. В то же время в современных инфокоммуникационных системах для 
решения функциональных задач используется большое число независимых генераторов [14—16]. При 
отсутствии высокостабильного генератора частоты каждого из совместно и независимо функциониру- 
ющих генераторов могут оцениваться так, как предложено в [17, 18]. Однако данный способ не позво- 
ляет учитывать изменения параметров генераторов при изменении внешних условий, что приводит к 
снижению точности оценок частот генераторов. 

Целью статьи является разработка метода, обеспечивающего повышение точности оценивания 
частоты генераторов в процессе жизненного цикла радиоэлектронной системы с учетом нестабильности 
температуры, напряжений и т. д. При этом предполагается, что отсутствует высокостабильный допол- 
нительный генератор и с этим связано непрогнозируемое отклонение длительности временного интер- 
вала. 

Теоретические основы метода оценивания частот генераторов в условиях непрогнозируе- 
мого изменения интервала измерений. Рассмотрим совокупность А+1 генераторов, объединен- 
ных в составе радиоэлектронной системы или комплекса. Номинальные значения частоты генераторов 





* Работа выполнена в рамках инициативной НИР. 
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и их относительные нестабильности известны и равны соответственно А, И а, (А=1,..,К +1). Для 
каждого из данной совокупности К генераторов в течение интервала длительностью 7, производится 
измерение полной фазы Ё колебаний его выходного сигнала ( К=1.....К’). Это происходит при 


текущих условиях эксплуатации, одновременно, в течение одного и того же временного интервала из- 
мерений, задаваемого (А’ + 1) -м генератором. 

В результате воздействия внешних факторов (например, нестабильности температуры, напря- 
жения), частоты всех генераторов отличаются от номинальных значений — на величину ДА, 
(К=1,...К +1), а длительность временного интервала — на величину ОТ, =-Т, х ЛЕ, / к... 

Чтобы оценить отклонение частоты каждого из совокупности А’ генераторов, измерим полную 
фазу его колебаний ф, и определим отклонение Дф, полной фазы колебаний относительно ее номи- 


нального значения ф,„ : 
ДАф‚ =Ф, — фи, К=1...К. (1) 

Значение ф., соответствует полной фазе колебаний сигнала номинальной частоты /), за вре- 
менной интервал измерений номинальной длительности 7,. 

Отклонение Лф, после линеаризации (отбрасывания слагаемого ДА, - 07, ) определяется двумя 
слагаемыми: 

— составляющей отклонения, определяемой отклонением частоты самого К -го генератора и 
равной Лф, = ДА -Т,; 

— составляющей отклонения, определяемой отклонением длительности временного интервала 
измерений от номинального значения и равной Дф, =№,-67,. 

При использовании метода ФАПЧ (К +1) -й генератор, применяемый для задания временного 
интервала измерений, является высокостабильным (его стабильность гораздо выше по сравнению с 
остальными А-генераторами). Поэтому Д/,., << ДА (К=1...., К), и для всех А-генераторов можно счи- 
тать Дф, << Дф, . С учетом последнего получаем, что Дф, = Дф, , и определяем отклонение частоты А 
го генератора ( А =1...., К) по результатам измеренного отклонения фазы колебаний от номинального 
значения с использованием соотношения ДА, = Дф, /Т, . Однако в рассматриваемом случае отсутствует 
высокостабильный генератор, обе составляющие отклонения фазы являются соизмеримыми и по ре- 
зультатам одиночных измерений не могут быть разделены. Таким образом, возникает задача оценки 
частоты генераторов в условиях непрогнозируемого изменения интервала измерений. 

Для разделения составляющих Лф, и Дф, необходимо получить оценку бт, и затем опреде- 
лить частоту колебаний каждого генератора при данных условиях эксплуатации. Для этого использу- 
ется соотношение: 

к_ к -д т 


ф 
Ё из 
^ 7 


и 


Для определения 67, составим функцию правдоподобия [19] 


Их (2) 


ит =] (по?) ^” ехр С, (3) 


2 
«а 20, 
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Коэффициенты /7, (К =1,..,К’) определяют степень влияния отклонения фазы каждого гене- 
ратора на точность оценивания отклонения длительности интервала измерений от номинального зна- 
чения. При записи соотношения (3) было учтено, что отклонения частот генераторов от номинальных 
значений ДА, (А =1,..., К’) определяются влиянием большого числа независимых и равнозначных фак- 


торов и, соответственно, подчиняются нормальному закону распределения. 
С учетом выражения (2) функция правдоподобия (3) приводится к виду 





К 1/2 ДАф, — 520 7 - 
бт} =т:] ] (2по?) " ехр ь ь / (4) 


В качестве искомой оценки 07, выберем значение, при котором функция (4) достигает макси- 
мума. Соответствующее значение равно 
К 
-2 р-1 
х У (т, Аф, ох ль ) 
= К=1 
ОТ, Т, г. : (5) 
Утв 
К=1 
Доказательство, что найденная В (5) оценка отклонения длительности временного интервала 


измерений обеспечивает максимум функции правдоподобия, легко получается из исследования второй 
производной функции (4). 





Значение 07, после подстановки в (2) позволяет получить оценку частоты А-го генератора 


(^=1,...К) с использованием выражения 


х Е. 
ети. (6) 


и 
Данное соотношение полностью решает поставленную задачу получения оценки частоты каж- 
дого из совокупности одновременно и независимо функционирующих генераторов в условиях непро- 
гнозируемого изменения интервала измерений. 
Введение коэффициентов т, позволяет увеличить число степеней свободы в управлении ве- 


ЛИЧИНОЙ получаемой оценки и компенсации воздействия различных факторов. В частности, как отме- 
чалось выше, в течение жизненного цикла происходит изменение параметров отдельных генераторов, 


изменяется их нестабильность 0”. Для компенсации указанных эффектов необходимо изменять пара- 
метры вычислительного алгоритма, на основе которых получается оценка 07, в (5). При этом доста- 


точно обеспечить выполнение Условия 1/7, ея = СОЛ5Ё. Это технически значительно проще, чем ис- 


пользование, например, систем термостабилизации, чтобы в течение всего цикла работы системы обес- 
печить стабильность частот генераторов. 
Свойства получаемых оценок. Исследуем свойства получаемых оценок. Для определения матема- 


^ 


тического ожидания получаемой оценки К, используем выражение 








Е | _ $ ть Е : б 7 
м ]=м| 9 (7) 
С учетом свойств математического ожидания преобразуем выражение (7) к виду 
1 М ]-^, -М[8Т, | 
м | Е - (8) 


и 
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Для дальнейших преобразований используем выражения (2) и (5). В частности, для М]ф, | 


можно записать 


Е 


М[$, |= фк +М[А*]-Т, ок ТМ Ай, 


Е 


0К-+1 
Из выражения (9) с учетом сделанных предположений о характере статистического распреде- 
ления отклонений частот генераторов непосредственно следует, что М [$ | = ф›. Аналогично можно 


(9) 





показать, что Л [67, | =0. Таким образом, из (8) следует, что математическое ожидание оценки ча- 


стоты А-го генератора равно №,, К=1....,К’, т. е. получаемая оценка частоты генератора является 
несмещенной [20]. 

Вычислим дисперсию получаемой оценки. Воспользуемся выражением (6), в котором случайной 
величиной является только оценка отклонения длительности временного интервала от номинального 
значения. В этом случае, основываясь на результатах [19], получим 


211-287.) К: (10) 


и 


Чтобы найти 215, | ‚ основываясь на (5) и учитывая, что первый начальный момент случайной 


величины 07, равен нулю, запишем выражение для второго центрального момента (дисперсии) оценки 
От, 
2 


ту, (о) 
267, = м Е (11) 


21>` 7,0 


К=1 





Выполняя последовательно преобразования в (10), получим 


г. кк кк р 
267, |=м $ у пытм,муорози 8; [5 у лытуаро? | (12) 


1 р=1 = р= 


Учитывая, что отклонения частот генераторов являются некоррелированными, получаем 


2{4$,} ‚ КЕр, 
м1, № |= 13 
[Аф, -ф,} В __ (13) 
Тогда представление для дисперсии бт, имеет вид 
я < 2 вх рН е 
р[6т,| = м (ть, (4) [УХ тьтьоро? |. (14) 
= К р— 


Принимая во внимание представление Дф, = Дф, + Дф, из (1), получаем 
2 
м (м, | = Во} +10). (15) 


Следовательно, окончательное выражение для дисперсии оценки временного интервала имеет 
вид: 


145 


Вестник ДГТУ. 2014. Т. 14, №4 (79) 








я —1 
кт (бе абка (9% 





216 7, = — (16) 
|} лутьовоз) 
К=1 р=1 
Таким образом, дисперсия оценки частоты генератора определяется формулой: 
ы 1 
Е о | 
и и (17) 
т [Хтьтьювоя) 
К=1 р=1 


Рассмотрим частный случай, при котором все генераторы работают на различных частотах 
Б.К (К=1....К +1) и имеют одинаковые относительные нестабильности а, =а (К =1....,К). Ис- 


следуем характеристики распределения отклонений частот генераторов от номинальных значений для 
данного случая. Исходя из того, что первый начальный момент случайной величины ДАф, равен нулю, 


на основании соотношения (16) запишем выражение для второго центрального момента (дисперсии): 
к 
2 +22 2 [2 2 
«О Ут, ПЕ ы пы) 
№ к 


кк 
и > тит, 


К=1 р= 


р (18) 





— 
=> 
ме 


В случае одинаковых значений весовых коэффициентов /л7, (К =1,....К) величина дисперсии 
определяется выражением 


К 
ое +) 
ы К? 





2{%} (19) 
Как следует из (18) 
т РА} =0 20 
К-—о { #} Е ( ) 
следовательно, данная оценка является асимптотически эффективной [20]. 
Для проверки свойства состоятельности получаемой оценки воспользуемся неравенством Чебы- 
шева. В соответствии с последним условием, для любого положительного числа = вероятность откло- 


нения оценки частоты генератора /, от его математического ожидания не меньше чем на = ограничена 
сверху величиной 2 {| /=?. 


В рассматриваемом случае получаем 





2] 

















Р(А >=) < к. (21) 
С учетом формулы (18) последнее выражение преобразуется следующим образом: 
с Вена. 
Р(А > =) < ТЕ : (22) 


откуда следует, что оценка является состоятельной [20]. 
Выводы. 

1. Предложенный метод статистической обработки результатов измерений фаз колебаний од- 
новременно и независимо функционирующих генераторов в условиях непрогнозируемого изменения 
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длительности интервала измерений позволяет получать оценки частот генераторов. При этом на пара- 
метры генераторов (частоту и стабильность частоты колебаний) не накладываются никакие ограниче- 
ния. Метод позволяет значительно упростить построение систем, обеспечивающих стабильность пара- 
метров генераторов, путем применения весовых коэффициентов, используемых при получении оценок 
частот генераторов. 

2. Как показали результаты исследований, получаемые с использованием данного метода 
оценки частот генераторов являются несмещенными, состоятельными и эффективными. 
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О. А. баГагуап 


Тре тейод о! п’едиепсу езйтайоп оЁ те хтийапеоиж/у апа таерепаепёйУу орегайпда депегаюг$ ипаег ипргеаа Ме 
свапдез п те теазигетепе /тёегиа! $е /5 виа ед. Еггог5 ап ипдег пе тетод арр/сайоп дие 0 те гапхепсу оЁ 
Те о5сюг Редиепсе$ оп Ве Редиепсу езИтайоп /тегиа! аге сопз/АЕгеа, ТПе В!5Ё еггог сотропепе /5 геаеа ю пе 
деиавоп оЁ пе теазига /е депегаог рпазе аие ю {те пвегепЕ депегаюг гедиепсу па; апа йе зесопа опе 15 
аегт/пед Бу пе ипргефдавае свапдез п {пе о5с//еог Редиепсу оп {пе теазигетепЕ /тегиа/. 1 15 поед паё те 
Десгеазе о! еасй о! йе сотропепёз /трозез псотраНЫЕе гедигетет оп пе дигайоп о {те йте !епдЁй. Оп ве Ва$5 
ОР {Пе кпоитп геайоп$ Гаепёйупд {пе роепёа/у геасва Ме иа/ие о! {пе гооЁ-теап-здиаге й’едиепсу аеиайоп от йе 
потГпа/ иа/ие, ве ехргеззоп5 ЧЕ егттиптд пе орйтит йте (епдЁ! аге обатеа, Аз а сйепоп иреп споо$пд те йте 
Гепа, те тигитит оЕ Вмо еггог$ 5ит /5 соп5/АЕгГед, Тре апа/уйса! ге/аНоп5 деептип!па роепва/у геаспа/е ассигасу 
ОР те бедиепсу езйтай=5 аге ргезетеа. Тре ипЫазеапес5, ейс/епсу, апа сопх&епсу оЁ 1Ше ета обеатеа аге 
ргоуед. 

Кеуигога5: депегаюг п’едиепсу, Ийте (еп? езйтайоп, еиамайоп оЁ Редиепсу демайоп тот пе пота! иа/ие, &а- 
В5Нса/ тетоа оЁ педиепс/е$ 5аб//ганоп, теазигетепеЕ /егиа/. 





* Тне гесеагсН {5 Чопе мил {Не #гате о! {Не (пдерепаепЕ ВО. 
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